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INTRODUCTION 


The objective of the present research is to develop a 
general mathematical model and solution methodologies for analyz- 
ing structural response of thin, metallic shell-type structures 
under large transient, cyclic or static thermomechanical loads. 
Among the system responses, which are associated with these load 
conditions, are thermal buckling, creep buckling and ratcheting . 
Thus, geometric as well as material— type nonlinearities (of high 
order) can be anticipated and must be considered in the develop- 
ment of the mathematical model. Furthermore, this must also be 
accommodated in the solution procedures. 


SUMMARY OF PROGRESS 

The progress to date has been elaborated upon in an interim 
scientific report submitted to the sponsor during the summer of 
1986, and in a series of semiannual progress reports. The most 
recent of these is dated April 1, 1987. 

A complete true ab-inito rate theory of kinematics and kine- 
tics for continuum and curved thin structures, without any res- 
triction on the magnitude of the strains or the deformation, was 
formulated. The time dependence and large strain behavior are 
incorporated through the introduction of the time rates of the 
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metric and curvature in two coordinate systems; a fixed (spatial) 
one and a convected (material) coordinate system- The relations 
between the time derivative and the covariant derivatives (gra- 
dients) have been developed for curved space and motion, so that 
the velocity components supply the connection between the equa- 
tions of motion and the time rate of change of the metric and 
curvature tenors. 

The metric tensor (time rate of change) in the convected ma- 
terial coordinate system is linearly decomposed into elastic and 
plastic parts. In this formulation, a yield function is assumed, 
which is dependent on the rate of change of stress, metric, tem- 
perature, and a set of internal variables. Moreover, a hypoelas- 
tic law was chosen to describe the thermoelastic part of the de- 
formation. 

A time and temperature dependent viscoplastic model was for- 
mulated in this convected material system to account for finite 
strains and rotations. The history and temperature dependence 
were incorporated through the introduction of internal variables. 
The choice of these variables, as well as their evolution, was 
motivated by phenomenological thermodynamic considerations. 

The nonisothermal elastic-viscoplastic deformation process 
was described completely by "thermodynamic state" equations. Most 
investigators (in the area of viscoplasticity) employ plastic 
strains as state variables. Our study shows that, in generalise 
of plastic strains as state variables may lead to inconsistencies 
with regard to thermodynamic considerations. Furthermore, the 
approach and formulation employed by all previous investigators 
lead to the condition that all plastic work is completely dis- 
sipated. This, however, is in contradiction with experimental 
evidence, from which it emerges that part of the plastic work is 
used for producing residual stresses in the lattice, which, when 
phenomenologically considered, causes hardening. Both limitations 
are not present in our formulation, because of the inclusion of 
the "thermodynamic state" equations. 

The obtained complete rate field equations consist of the 
principles of the rate of the virtual power and the rate of 
conservation of energy, of the constitutive relations, and of 
boundary and initial conditions. These formulations provide a 
sound basis for the formulation of the adopted finite element 
solution procedures. 

The derived shell theory, in the least restricted form, 
before any simplifying assumptions are imposed, has the following 
desirable features: 

(a) The two-dimensional, impulse-integral form of the equations 
of motion and the Second Law of Thermodynamics (Clausius- 
Duhem inequality) for a shell follow naturally and exactly 
from their three-dimensional counterparts. 

(b) Unique and concrete definitions of shell variables such as 
stress resultants and couples, rate of deformation, spin and 


entropy resultants can be obtained in terms of weighted 
integrals of the three-dimensional quantities through the 
thickness . 

(c) There are no series expansions in the thickness direction. 

(d) There is no need for making use of the Kirchhoff Hypotheses 
in the kinematics. 

(e) All approximations can be postponed until the descretization 
process of the integral forms of the First Law of Ther- 
modynamics 

(f) A by-product of the descent from three-dimensional theory 
is that the two-dimensional temperature field (that emerges) 
is not a through-the-thickness average, but a true point by 
point distribution. This is contrary to what one finds in 
the literature concerning thermal stresses in the shell. 

To develop geometrically nonlinear, doubly curved finite 
shell elements the basic equations of nonlinear shell theories 
have to be transferred into the finite element model. As these 
equations in general are written in tensor notation, their 
implementation into the finite element matrix formulation 
requires considerable effort. 

The nonlinear element matrices are directly derived from the 
incrementally formulated nonlinear shell equations, by using a 
tensor-oriented procedure. The classical thin shell theory based 
on the Kirchof f-Love hypotheses (Formulation D) was employed for 
this purpose. For this formulation, we are using five "natural" 
degrees of freedom per mid-surface shell node: three incremental 
velocities and two rates of rotations about the material coor- 
dinates «* 1 , and<* 2 , in a mixed form. 

The quasi-linear nature of the principle of the rate of 
virtual power suggests the adoption of an incremental approach to 
numerical integration with respect to time. The availability of 
the field formulation provides assurance of the completeness of 
the incremental equations and allows the use of any convenient 
procedure for spatial integration over the domain V. In the 
present instance, the choice has been made in favor of a simple 
first order expansion in time for the construction of incremental 
solutions from the results of finite element spatial integration 
of the governing equations. 

The procedure employed permits the rates of the field formu- 
lation to be interpreted as increments in the numerical solution. 
This is particularly convenient for the construction of incremen- 
tal boundary condition histories. 

Even under the condition of static external loads and slowly 
growing creep effects, the presence of snap-through buckling 
makes the inertia effects significant. In dynamic analyses, the 
applied body forces include inertia forces. Assuming that the 
mass of the body considered is preserved, the mass matrix can be 
evaluated prior to the time integration using the initial 
configuration. 


Finite element solution of any boundary-value problem 
involves the solution of the equilibrium equations (global) 
together with the constitutive equations (local) . Both sets of 
equations are solved simultaneously in a step by step manner. The 
incremental form of the global and local equations can be 
achieved by taking the integration over the incremental time step 
&t=tj + 1 -tj . The rectangular rule has been applied to execute the 
resulting time integration. 

Clearly, the numerical solution involves iteration. A simp- 
lified version of the Riks-Wempner constant-arc-length method has 
been utilized. This iteration procedure which is a generalization 
of the displacement control method also allows to trace the non- 
linear response beyond bifurcation points. In contrast to the 
conventional Newton— Raphson techniques, the iteration of the 
method takes place in the velocity and load rate space. The load 
step of the first solution in each increment is limited by 
controlling the length ds of the tangent. Either the length is 
kept constant in each step or it is adapted to the characteris- 
tics of the solution. In each step the triangular-size stiffness 
matrix has to be checked for negative diagonal terms, indicating 
that a critical point is reached. 

One of the most challenging aspects of finite strain formu- 
lations is to locate analytical solution with which to compare 
the proposed formulation. Typically, as a first problem, a large 
strain uniaxial test case was analyzed. The case considered 
examines the rate-dependent plastic response of a bar to a 
deformation history that includes segments of loading, unloading, 
and reloading, each occurring at varying strain and temperature 
rates. Moreover, it was shown that the proposed formulation 
generates no strain energy under a pure rigid body rotation. 
These are surely important demonstration but they only represent 
a partial test because the principal stretch directions remain 
constant. Finally, a problem which was discussed by Nagtegaal and 
de Jong, and others too, as a problem which demonstrates limita- 
tions of the constitutive models in many strain formulation, is 
the Couette flow problem. This problem is solved as a third 
example. The results of these test problems show that: 

The formulation can accommodate very large strains and 
rotations . 

The formulation does not display the oscillatory beha- 
vior in the stresses of the Couette flow problem. 

— The model incorporates the simplification associated 

with rate-insensitive elastic response without losing 
the ability to model rate temperature dependent yield 
strength and plasticity. 

The problem of buckling of shallow arches under transient 
thermomechanical load was investigated next. The analysis was 


performed with the aid of 24 paralinear isoparametric elements. 
The paralinear isoparametric element is such that the thickness 
is small compared to other dimensions. The characteristics of the 
element are defined by the geometry and interpolation functions, 
which are linear in the thickness direction and parabolic in the 
longitudinal direction. Consequently, the element allows for 
shear strain energy since normals to a mid-surface before 
deformation remain straight, but not necessarily normal to the 
midsurface after deformation. 

The developed solution scheme is capable of predicting 
response which includes pre- and post-buckling with thermal creep 
and plastic effects. The solution procedure was demonstrated 
through several examples which include both creep and snap- 
through behavior. 

The last set of problems which are under investigation 
consists of creep or thermal buckling, with plastic effects, of 
shells of revolution. 

In addition, following a more traditional approach, a method 
was developed for bounding the response (solution) of bars and 
beams of (linear) viscoelastic material behavior, based on 
nonlinear kinematic relations. 

In connection with this progress to date, two papers were 
published by the AIAA Journal in 1986. Moreover, a paper en- 
titled, "Non-Isothermal Elastoviscoplastic Analysis of planar 
Curved Beams" was presented at the 3rd Symposium on Nonlinear 
Constitutive Relations for High-Temperature Applications, held at 
the University of Akron, on June 11-13, 1986. A descriptive 
abstract of this paper was published in the meeting proceedings 
and the full paper will appear in a special publication. Copies 
of the above have been sent to the sponsor. 

In addition, two papers were presented at the 28th AIAA/- 
ASME/ASCE/AHS SDM Conference and published in the proceedings of 
this conference, and one paper was presented at the 21st Con- 
ference of the Israel Society of Mechanical Engineers (keynote 
lecture by Dr. R. Riff) . All three papers have been accepted for 
publication, the first two by the AIAA Journal and the last one 
by the International Journal of Computers and Structures. These 
three papers deal with applications to snap-through and creep 
buckling of bars and arches, and to two-dimensional problems in 
extension and shear. Most of this work was also presented at the 
NASA-Lewis Conference on Structural Integrity and Durability of 
Reusable Space Propulsion Systems on May 1987 in Cleveland. 

Progress has also been made with applications to shell 
configurations, and papers are in preparation dealing with these 
systems . 


FUTURE TASKS 


At this stage, the study is limited only to the simplest of 
the developed shell theory formulation (Formulation D) . 

The main thrust of the additional tasks is to develop a 
finite element and select a code, which will be made available to 
all interested users and which will be based on the most general 
nonlinear shell formulation possible and nonlinear constitutive 
relations to predict the response of shell-like structures, when 
subjected to time-dependent thermomechanical loads with large 
excursions. 

Before outlining the research program, objectives and tasks, 
a brief description of the various nonlinear shell theories is 
presented next (for more details see the 1986 technical report) . 

The various shell theory approximations (formulations) are 
based on the use of certain simplifying assumptions regarding the 
geometry and kinematics of the shell configuration. 

These are: 

Assumption I : The material points which are on the normal to 
the reference surface before deformation will be on the same 
normal after deformation. 

Assumption II : The shell is sufficiently thin so that we can 
assume linear dependence of the position of any material 
point (in the deformed state) to the normal (to the refe- 
rence surface) coordinate (in the deformed state) . The 
linear dependence can easily be changed to parabolic, cubic, 
or any desired degree of approximation. 

Assumption III : The rate of change of the velocity gradients 
with respect to in-plan coordinates on the two bounding 
shell surfaces is negligibly small. 

Assumption IV : The rate of change of the distance of a 

material from the reference surface is negligibly small. 

On the basis of the above four simplifying assumptions, 
several formulations result, for the analysis of thin shells. 
These formulations are denoted below by capital letters. 

Formulation A : This formulation makes use of Assumption 

I, only . 

Formulation B : This formulation employs Assumptions I and 

II. 

Formulation C : This formulation employs Assumption I, II and 

IH. 

Formulation D : This is the classical thin shell theory based 
on the Kirchhoff-Love hypotheses of Assump- 
tions I, II, III, IV, as applied to large 
deformation theory. 

These formulation are arranged in such a manner that we move from 


the least restrictive (A) to the most restrictive (D) . 


In addition to this a fifth formulation (E) can easily be 
devised and this formulation in terms of order of restriction is 
similar to Formulation A. Formulation E makes use of Assumption 
II only. 

In order to accomplish the overall objective, as stated ini- 
tially, a number of steps or tasks need be completed. These are: 


Step 1 : 


Step 2 : 
Step 3 : 


Step 4 : 


Step 5 : 
Step 6 : 


Formulation the principle of the rate of virtual 
power for the less restricted shell formulation 
(C, B, E, A). 

Couple the general constitutive relations to shell 
kinematics for the aforementioned formulations. 
Before developing a finite element, and before 
numerically testing and comparing the merits of 
the various formulations, some theoretical 
evaluation of the degree of approximation will be 
attempted. 

Develop a finite element for the chosen least 
restrictive shell formulation. In order to 
accomplish this some novel ideas and procedures 
will be employed. These include a tensor-oriented 
procedure for obtaining the element, and a 
discretization process that involves not only the 
nodal degrees of freedom but also an approximation 
to the "shift” tensors. 

Incorporate the developed finite element and 
related constitutive relations into a FEM code. 
Numerically test and evaluate the developed finite 
element procedure. 


It is estimated that this entire research effort will 
require a minimum of three years to bring to fruition. Some of 
the tasks (steps 1 and 2) have been started and will be continued 
during the current funding year. 


